Microscopic origin of Magnetic Ferroelectrics in Nonlinear Multiferroics 
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A simple but general microscopic mechanism to understand the interplay between the electric and 
magnetic degrees of freedom is developed. Within this mechanism, the magnetic structure generates 
an electric current which induce an counterbalance electric current from the spin orbital coupling. 
When the magnetic structure is described by a single order parameter, the electric polarization is 
determined by the single spin orbital coupling parameter, and the material is predicted to be a half 
insulator. This mechanism provides a simple estimation of the value of ferroelectricity and sets a 
physical limitation as well. 



Multiferroics are materials in which magnetic and elec- 
tric orders are strongly coupled, and have attracted in- 
creasing attention^'^. Recent experimental research on 
multiferroics has shown that ferroelectricity and mag- 
netism not only coexist in the same material but also 
couple so strongly that the magnetic degree of freedom 
can be manipulated by an electric field and the elec- 
tric degree of freedom can be manipulated by a mag- 
netic field^i^'^i^'^'^i^. This property promises important 
technological applications in the future. However, the 
strong coupling between these two degrees of freedom 
poses enormous challenges to theoretical attempts to un- 
derstand the microscopic mechanism inside the materi- 
als. So far, due to the complexity of such materials, the- 
oretical understanding of the microscopic mechanism of 
strong magnetoelectric coupling is quite limited. 

Through symmetric analysis in the Ginzburg-Landau 
approach to thermodynamics, M. Mostovoy^° has shown 
that, the relationship between the ferroelectric order P 
and magnetic order M in spiral magnets is given by 



P cx Af X (V X M). 



(1) 



Eq.l qualitatively explains experimental results. Mi- 
croscopically, a mechanism, an inverse effect of the 
Dzyaloshinskii-Moriya (DM) ^^'^^ interaction, to generate 
EqUhas been proposedi^. The inverse DM mechanism is 
based on the idea that spin currents are induced between 
the noncoUinear spins, and can therefore be considered 
as electric moments. This is the only known mechanism 
to date. However, this mechanism is limited in two im- 
portant ways: (1) it lacks quantitative prediction results; 
(2) strictly speaking, the concept of spin current in spin 
orbital coupled system is not well defined. This Letter 
proposes a simple but general new mechanism to under- 
stand Eq.l. Based on well defined conventional concept, 
electric current, the proposed mechanism is completely 
different from the inverse DM mechanism. The proposed 
mechanism not only can provide quantitative prediction 
but also enable us to consider a unified picture of strong 
magnetoelectric coupling in multiferroics. 

We know from Maxwell equation that static magnetism 
and static electricity do not couple with each other, and 
that dynamic electricity such as an electric current can 
generate a magnetization. The inverse process also exists. 



namely, a magnetization can generate an electric current. 
Let's revisit the definition of the electric current. The 
current operator of the electron is defined as the change 
in Hamiltonian with respect to the variation of the vector 
potential of electromagnetic field, i.e. 
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(2) 



In non-relativistic quantum mechanics, the definition of 
the electric current includes three terms generated from 
three different physics: (1) the contribution of standard 
momentum; (2) the spin contribution illustrated in stan- 
dard quantum textbooks-^^; (3) the contribution of spin 
orbital coupling. To be more specific, we consider a single 
electron in a band structure. The electron is described 
by the Hamiltonian 



A\2 



2m* 



a{p — e — ) 
c 



(ct X VV{f)) - X A) ■ ct(3) 



the 



where, m* is the effective mass of electrons, a is 
effective spin orbital coupling parameter, /x = 2^ and a 
is the spin of the electron. In the absence of the external 
electrodynamic field, i.e. A = 0, for a given wavefunction, 
\l/(r), the electric current derived from Eql2]is given by 

j ^ Jo + A'cV X + Q;e(^'*CT^') X Vy(f) (4) 



where 



Jo 



2m' 



(5) 



The three terms in Eq.4 precisely correspond to the three 
physics discussed above. The first term usually domi- 
nates over the other two terms in transport. Therefore, 
the last two terms are usually ignored and are not famil- 
iar to most people. Here we show that in multiferroics 
which are insulators, the interplay between the last two 
terms in Eqd] provides a fundamental mechanism to un- 
derstand the magnetoelectric coupling. 

In the nonlinear multiferroics, RMnO^, a nonlinear 
spiral magnetic order has been observed^. The spiral 
magnetic order is formed by the localized spins of Mn 
atoms. To study the electronic physics, we can use a 
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Kondo-lattice type of model. We consider that the elec- 
trons in the band couple to the localized spins of the Mn 
atoms through spin exchange coupling. Through the ex- 
change coupling, we can naturally assume that the mag- 
netic ordering of the localized spins also generates the 
same magnetic ordering for the electrons in the band. 
Now consider that the magnetization of the electrons in 
the band is a simple spiral magnetic ordering 



Mq = M{){cosqx/a, sinqx/a, 0). 



(6) 



The electric current associated with the magnetization is 
given by 



Jm = AicV X Mo — 



^cqMo 



(0, 0, cosqx/a) 



(7) 



which is a current along z direction. In fact, it is a 'global' 
current along z direction for a fixed x coordinate. In an 
insulator, the net electric current with such a configura- 
tion must be zero based on Kohn's proof of the insulat- 
ing property namely, the magnetization current must 
be counterbalanced by other electric currents. The total 
electric current contributed from jo in the band also van- 
ishes since the lattice mirror symmetry in the x-y plane 
is not broken in the nonlinear multiferroics^ii^iii in the 
absence of external magnetic field. Therefore, the elec- 
tric current from the magnetic ordering must be coun- 
terbalanced by the electric current induced from the spin 
orbital coupling. 

The above analysis can also apply to general nonlinear 
magnetic ordering structures which induce similar elec- 
tric currents. From Eq|31 the cancelation requirement 
leads to 



/xcV X Mo + aeMo x W{r) = 



(8) 



By simple algebraic modification and averaging in the 
total space, Eq[5]beomes 



< E >=< 



(Mo • E)Mo 



Mi 



> + < ^^^^^^ >(9) 



where < ... > takes the space average and W{r) = 
—eE{r). The first term in the right side of above equa- 
tion usually vanishes when taking the space average for 
a space modulating spin density. We obtain the total 
ferroelectricity as 



- _ eoMC Mq X (V X Mo) 
ae2 ^ M2 



(10) 



EqfTOl is consistent with EqlTJ However, Eq[TO] pro- 
vides detailed coupling coefRcients that are different from 
the results normally expected from the Ginzburg-Landau 
theory. First, the coefficients are inversely proportional 
to effective spin orbit coupling, which is against intuitive 
expectation. Second, Eg 11 01 predicts that saturated value 
of ferroelectricity does not depend on the amplitude of 
the magnetization. Therefore, it suggests that increasing 



large magnetization will not dramatically increase ferro- 
electricity, which is against that normally expected from 
the simple Ginzburg-Landau Theory. Finally, in EqfTUl 
there is only one free parameter, the effective spin or- 
bital coupling. The fact of the existence of only one free 
parameter makes relatively easier to test the new mech- 
anism experimentally. 

There are two important issues regarding of the above 
results. First, it is clear that the above counter-intuitive 
dependence of the polarization on the spin-orbit coupling 
parameter and the magnitude of magnetization can not 
be correct for arbitrary small values of these parameters. 
Therefore, what is the limitation of this mechanism based 
on the current cancelation? To answer this question, we 
have to compare the energy saved from the current can- 
celation and the energy cost from the polarization. When 
a polarization is developed, there is an energy paid for 
the deformation of the lattice. This energy cost Ecost for 
a small polarization is expected to be 



(11) 



where A depends on the detailed lattice structures. The 
energy saved in the current cancelation for the electrons, 
Esave can also easily estimated. 



(12) 



In order to favor the current cancelation mechanism, we 
must have Esave > Ecost, which leads to the following 
criteria for the mechanism of the current cancelation by 
plugging EgfTOl 



aMo > 



(13) 



This criteria is satisfied in the spiral magnets such as 
TbMnOs as we will show later. 

Second, a careful reader may notice that there is a 
critical flaw in the above derivation for a simple band 
picture: for any magnetic ordering of the localized spins 
of atoms, the electrons in a fully filled band do not have 
magnetization response regardless of the strength of spin 
exchange coupling between them. Namely, Mo = for 
any magnetic order M of the localized spins of atoms. 
Thus the cancelation of the current in Eq[5] does not ex- 
ist. This observation leads to another important predic- 
tion of the paper: in order to generate magnetoelectric 
coupling through the current cancelation, the multifer- 
roics must be a 'half insulator. Let's consider a single 
spin degenerate band to understand the rational behind 
the prediction. The presence of electron-electron interac- 
tion and the exchange interaction between the electrons 
in the band and the localized spin moments can all cause 
a single spin degenerate band to split into several bands. 
The simple picture in Fig.l shows the splitting of the 
original bands in two. If both of the two new bands 
are fully occupied, no electric current will be generated 
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by the magnetic ordering of the locaHzed spins, because 
Mq = 0. As a matter of fact, one can picture the physics 
as the magnetization and electric currents are exactly the 
opposite in direction in two new bands. However, if a gap 
exists between the two new bands, and if the lower en- 
ergy band is completely filled but the upper energy band 
is empty, which we call it as a 'half insulator. The curl 
of the magnetic ordering generates real electric currents 
in the half insulator which needs to be counterbalanced 
by the electric current from the spin orbital coupling. 

Now we quantitatively discuss the ferroelectricity pre- 
dicted by Eqlini We rewrite a = where can 
be viewed as the effective spin split energy gap. For the 
spiral magnetic structure in EqlBJ if we choose reasonable 
parameters with g factor to be two, the lattice constant 
a = Inm and As with an unit of ev, we obtain 

Py ^ 0.88qAs[fiC/{cm^ ■ ev)] (14) 

If we make reasonable assumption that As varies from 
O.lev to lev and consider the fact that theoretical values 
are generally larger than experimental measurements due 
to the existence of disorder, we can conclude that values 
predicted by Eg 1141 are in the same order of experimental 
values observed in experiment o^^i^^ for the spiral magnets 
such as TbMnOs and DyMnOz. 

To show that the condition, Eq[T3l is really satisfied in 
the spiral magnets, we can estimate the cost energy Ecost 
and the saving energy Esave- Taking the polarization 
value Py ^ Q.\^C/cm? observed in TbMnO^, we can 
estimate the effective lattice shift, Sy = PyVn/e*, where 
Vn is the volume of one unit cell and e* = N^ffe is the 
effective charge in one unit cell. The cost energy Ecost = 
^KSy where k is determined by lattice structures. For 

TbMnOs, K - lev/A^ and K = 2.3 x 10^1^, we have 
Ecost ~ ( ^°^^ )^(ei;). Let's assume that the number of 
effective spin of electrons in one unit cell Mq = Seff{fi), 
m* — bm (m electron mass) and g—2, we obtain, Egave ^ 
Q.ibbSlff lev) for g - 27r x 0.3 in TbMnOs. Comparing 
Ecost and Egave, we see that the EqUJlcan be satisfied for 
Scff > 0.02. It is important to note that the Seff is not 
the total effective spin of magnetized atoms measured in 
experiments and it is always less than one since there is 
no contribution from completely filled bands. 

Although the condition, Eq. [131 posts a limitation on 
the current cancelation mechanism, Eq [HI suggests that 
the value of ferroelectricity can still grow largely as As 
increases. However, in the electric current cancelation 
mechanism, there is an additional limitation on the value 
of the polarization. The ferroelectricity is limited by a 
natural energy scale, the energy gap Ag in the 'half in- 
sulator. The gap defines a length scale I — , ^ In 

order to maintain the validity of a 'half insulator, the 
value of ferroelectricity must satisfy the following crite- 



ria, 

|P| < =0.72,/6A3[/iC/(cm2.ei;i)] (15) 

el V " 

the energy unit of Ag is ev. Eq[T5] shows that it is 

very hard for the ferroelectricity in the multiferroics to 

be larger than a few fjiC/cm? because Ag should be 

in a range of a few ev at most. This explains why 

the ferroelectricity measured in the multiferroics is much 

lower than that in conventional ferroelectric materials. 

For example, the largest ferroelectricity in the family of 

RMnOa, R = Dy, Tb, Gd... is P ~ O.i^iC / cm'^^'^ which is 

measured in DyMnO^ while in conventional ferroelectric 

materials, such as BaTiOs, P = 26fiC/cm'^--. 

In summary, this letter develops a new mechanism to 
explain magnetoelectric coupling. With a relative small 
number of free parameters, the mechanism can estimate 
the value of ferroelectricity. The free parameters can be 
independently measured using various experimental tech- 
niques. For example, the energy gap in the multiferroics 
and the spin orbit coupling strength can be measured in 
optical absorption spectrum. The prediction of a 'half 
insulator can be tested via numerical calculations as well. 
Most importantly, the mechanism sets a general guide- 
line for the search of new multiferroics materials with 
ferroelectricity of larger value. 

The new mechanism also provides a new perspective 
through which magnetoelectric coupling can be under- 
stood: focusing on electronic physics rather than on lat- 
tice dynamics. The lattice dynamics has been the key 
to understand conventional ferroelectricity materials be- 
cause the development of ferroelectricity is tied to phonon 
softeningi^. However, in many multiferroics, no clear 
indication for the phonon softening has been found^i^^. 
Our new mechanism clearly states that in the multifer- 
roics, it is the electronic properties that are responsible 
for the magnetoelectric coupling. 

It is natural to ask if the mechanism we discussed above 
is universal to all multiferroics. Although it is not easy to 
test, we believe that the answer is positive. The analysis 
in this letter assume the presence of a single magnetic 
order in the material. With a single magnetic order, the 
electric current associated to the magnetic order is easy 
to be defined. In many multiferroics, the magnetic struc- 
ture are very complicated and can not be described by a 
single magnetic order parameter— i^. A careful construc- 
tion of electric current due to magnetization is required. 
These physics await future investigations. 
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FIG. 1: A sketch of two bands which are spht from one spin 
degenerate band. Ag defines the band gap. When the chem- 
ical potential at Ei, both bands are fully filled. When the 
chemical potential at E2, it is a 'half insulator which can 
generate ferroelectricity. 
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